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Abstract. We show that if the dwarf-nova disc instabil- 
ity model includes the effects of heating by stream impact 
and tidal torque dissipation in the outer disc, the calcu- 
lated properties of dwarf-nova outbursts change consider- 
ably, and several notorious deficiencies of this model are 
repaired. In particular: (1) outside- in outbursts occur for 
mass transfer rates lower than in the standard model as 
required by observations; (2) the presence of long (wide) 
and short (narrow) outbursts with similar peak luminosi- 
ties is a natural property of the model. Mass-transfer fluc- 
tuations by factors 2 can explain the occurrence of both 
long and short outbursts above the cataclysmic variable 
period gap, whereas below 2 hr only short normal out- 
bursts are expected (in addition to superoutbursts which 
are not dealt with in this article). With additional heating 
by the stream and tidal torques, such fluctuations can also 
explain the occurrence of both outside-in and inside-out 
outbursts in SS Cyg and similar systems. The occurrence 
of outside-in outbursts in short orbital-period, low mass- 
transfer-rate systems requires the disc to be much smaller 
than the tidal-truncation radius. In this case the recur- 
rence time of both inside-out and outside-in outbursts 
have a similar dependence on the mass-transfer rate M2. 
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1. Introduction 

Dwarf novae (DN) are erupting cataclysmic variable stars. 
The 4-6 mag outbursts last from few days to more than 
a month and the recurrence times can be as short as a few 

In 



days and as long as 30 years (see e.g. [Warner 1995[ ). 
cataclysmic variables, a Roche-lobe filling low-mass sec- 
ondary star is losing mass which is accreted by a white 
dwarf primary. In DNs the matter transfered from the sec- 
ondary forms an accretion disc around the white dwarf; 
it is this disc which is the site of the outbursts. It is 



believed that an accretion disc instability due to par- 
tial hydrogen ionization is triggering the outbursts. The 
disc instability model (DIM; see Cannizzo 1993, Lasota 
2000 for reviews) is supposed to describe the whole DN 
outburst cycle. This model assumes the a-prescription 



(Shakura & Sunyaev 1973) for the viscous heating and the 
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angular momentum transport in the disc. In such a frame- 
work, the outbursts are due to propagating heating and 
cooling fronts, while during quiescence the disc is refill- 
ing the mass lost during the outburst. In the standard 
version of the model, the mass transfer rate is assumed 
to be constant during the whole outburst cycle. In order 
to reproduce outburst amplitudes and durations, the a- 
parameter has to be larger in outburst than in quiescence 
(Smak 1984). 

The DIM is widely considered to be the model of DN 
outbursts because it identifies the physical mechanism giv- 
ing rise to outbursts, and with the above mentioned as- 
sumptions it produces outburst cycles which roughly look 
like the real thing. A closer look, however, shows that the 
model fails to reproduce some of the fundamental proper- 
ties of dwarf-nova outbursts. The main deficiencies of the 
DIM were recently summarized by Smak (2000). First, 
models predict a significant increase of the disc luminos- 
ity during quiescence (see e.g. Fig. an effect which is 
not observed. Second, in many systems the distribution of 
burst widths is bimodal (van Paradijs 1983). The DIM can 
reproduce the width bimodality but, contrary to observa- 
tions, its narrow outbursts have lower amplitudes than 
the wide ones. Third, observations clearly show that dur- 
ing outbursts the secondary is irradiated and the mass 
transfer rate increases. The standard DIM does not take 
this into account. To these problems pointed out by Smak 
one can add several other deficiencies. The model predicts 
outbursts beginning in the outer disc regions ('outside-in') 
only for high mass transfer while observations show the 
presence of such outbursts at apparently low mass transfer 
rates. In most cases the model predicts quiescent accretion 
rate which are at least two orders of magnitude lower than 
the rates deduced from quiescent luminosities. Finally, it 
is clear that superoutbursts observed in SU UMa stars 
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and other systems (Smak 2000) cannot be explained by 
the standard DIM. 

All this leads to the obvious conclusion that the stan- 
dard version of the DIM has to be enriched by inclu- 
sion of physical processes which might be important but 
which, for various reasons, were not taken into account in 
the original model. Osaki (1989) proposed that superout- 
bursts are due to a tidal-thermal instability. Smak (2000) 
suggested that problems of this model could be solved by a 
hybrid which would combine it with the radiation-induced 
enhanced mass transfer model (Osaki 1985; Hameury, La- 
sota & Hure 1997). He included mass transfer enhance- 
ments during outburst into the standard model and sug- 
gested that this could explain properties of the narrow and 
wide outbursts which alternate in many systems (Smak 
1999a). Hameury, Lasota & Warner (2000), extending pre- 
vious work by Hameury et al. (1997) and Hameury, Lasota 
& Dubus (1999), showed that inclusion of effects such as 
mass transfer variations, formation of holes in the inner 
disc and irradiation of both the disc and the secondary 
allows to reproduce many properties of the observed out- 
bursts. 

Some questions, however, were left without satisfac- 
tory answers. First, the problem of 'outsidc-in' outbursts 
occurring for too high mass-transfer rates was still un- 
solved. Second, Smak (1999a) obtained narrow and wide 
outbursts by including irradiation induced mass trans- 
fer fluctuations but wide outbursts in his model result 
from fluctuations whose amplitude bring the disc to a 
stable state. In such a case the width of wide outbursts 
would not be determined by the disc structure alone. 
There is nothing fundamentally wrong with such a model 
and long standstills of Z Cam systems are most prob- 
ably due to such fluctuations, as proposed by Meyer & 
Meyer-Hofmeister (1983) and shown by Buat-Menard et 
al. (2000). However, one should check that the model is 
complete, that no other mechanism can explain outburst 
width distribution. 

In the present paper we study two mechanisms which 
dissipate energy in the outer parts of an accretion disc. 
They are: (1) the impact on the outer boundary of 
the accretion disc of the gas stream which leaves the 
secondary through the Li critical point of the Roche 
lobe; this impact zone is observed in the light-curves 
as a bright ('hot') spot; (2) the viscous dissipation in- 
duced by the tidal torques which are responsible for 
the disc outer truncation. These mechanisms which may 
heat up the disc's outer regions have been known for 



a long time ( 


Paczyhski 1977 




Papaloizou & Pringle 1977 


Meyer & Meyer-Hofmeistcr 1983 


), but have been ne- 



glected in most of DIM calculations. The exception is 
Ichikawa & Osaki (1992) who included both effects in their 
model. They used, however, a viscosity prescription spe- 
cially designed to suppress inside-out outbursts and did 
not perform a systematic study of the effects as a func- 
tion of the mass-transfer rate, size of the disc etc. 



We add stream-impact and tidal torque heating effects 
to the Hameury et al. (1998) version of the DIM, pre- 
sented in Sect. 2. We use a ID approximation because the 
dynamical time-scale is shorter than the thermal time- 
scale. We study the effect of the stream impact (Sect. 3) 
and the tidal torques (Sect. 4) on outburst properties in 
short and long period DN systems. Both effects change 
the energy equation of the standard DIM. In both cases 
we find similar results: when compared to the standard 
case, 'outside-in' outbursts occur for lower mass transfer 
rate. The recurrence time is decreased by both the stream 
impact and dissipation of the tidal torques. We conclude, 
therefore, that heating of the outer disc by the stream 
impact and the rate of work by tidal torques provides 
a satisfactory explanation of the presence of 'outside-in' 
outbursts at realistic mass-transfer rates. 

The most important result of our investigation is that 
when these effects are included one can obtain two types 
of outbursts: short ('narrow') and long ('wide') (Sect. 5). 
We find a critical mass transfer rate for which these two 
type of outburst 'alternate'. This mass transfer rate is very 
close to that deduced from evolutionary calculations for 
systems with orbital periods longer than 3 hours ('above 
the period gap'). For lower mass-transfer rates only nar- 
row outbursts are present; for higher ones, outbursts are 
only of the wide type. Fluctuations of the mass-transfer 
rate allow to expand the range of orbital periods for which 
the bimodal outburst width distribution is occurring. We 
conclude that the observed bimodal outburst width dis- 
tribution is caused by the outer-disc heating. 

In Sect. 6 we discuss conditions for which outside-in 
outbursts can occur. We show that one can naturally ac- 
count for the occurrence of both types of outbursts in SS 
Cyg if the mass transfer rate fluctuates by factors < 2 
around a mean value which is close to what is generally 
assumed as the average mass transfer rate in this system. 
Outside-in outbursts in systems below the period gap, 
which have low mass transfer rates, can be explained only 
if the disc is significantly smaller than the tidal truncation 
radius. We show that this is in agreement with observa- 
tions and consistent with expectations of the tidal-thermal 
instability model, at least during the first half of the super- 
cycle of SU UMa systems. We also point out the ambiguity 
in observational determination of outburst 'outside-in' or 
'inside-out' type, also respectively called A and B types. 
These types are generally assigned to outbursts according 
to their shape and the UV delay. We question the general 
validity of such criteria and argue that the eclipse profile 
observations during outbursts is the only reliable diagnos- 
tics of the outburst type. 

2. The model 

We use the recent version of the disc instability model 
described in Hameury et al. (1998). We assume here 
that a = ttcoid = 0.04 in quiescence and a = ahot = 
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0.2 in outburst ( ^mak 1999b ). Difficulties in reproduc- 
ing the observed UV delays (see Sect. 6 for more de- 
tails) and accounting for observed quiescent X-ray fluxes 
well above the values predicted by the DIM, led to the 
suggestion that 'holes' are present in the center of qui- 



Table 1. Parameters adopted for long and short period 
systems. 



escent accretion discs (Meyer & Meyer- Hofmcister 1994, 
Lasota, Hameury fc Hurc 1995| ). Smak (1998) pointed 
out, however, that the failure to reproduce observed UV 
delays by numerical codes was largely due to the as- 
sumption of a fixed outer disc radius, which prevented 
the occurrence of outside-in outbursts. In addition, the 
discs used in such calculations were too small to give 
the required time-scales. When one uses reasonable outer 
boundary conditions, long UV delays can be observed. 
However, if one wishes to reproduce both the delays and 
the light-curve shapes, a truncated disc might be still a 
necessary ingredient of the model (Hameury et al. 1999). 
Such a hole could be due to magnetic disruption in the 
case of a magnetic white dwarf or to evaporation in the 
general case. Because we are not interested here in the 
exact values of the rise time and the UV delay, we use a 
fixed size hole only for numerical convenience. The inner 
disc radius is fixed at Tin = 10^ cm, while 5 x 10* cm is 
approximately the radius of a 1.2 M© white dwarf and 
6.9 X 10® cm correspond to the other mass we use: 0.8 
Mq . The outer boundary condition depends on the circu- 
larization radius Tcirc at which a particle that leaves the 
Lagrangian point Li with angular momentum j — fJorb^^ 
(where ilorb is the orbital angular velocity and b the dis- 
tance between the primary and the Li point), would stay 
in circular orbit if there were no accretion disc. From 
j = r^K^'circ' where I^k is the Keplerian angular frequency, 
one gets: rdrc — ^orb^''/^^^- When one takes the grav- 
itational potential of the secondary into account, these 
relations are slightly modified, and we linearly interpolate 
the table from Lubow & Shu (1975) to calculate rcirc- The 
mean outer radius < Tout > is taken to be the average 
of the three radii ri , r2 and r^^ax calculated in table 1 of 
Paczyirski (1977). 

We consider in this paper the effects of the impact of 
the gas stream coming from the secondary (Sect. 3) and of 
the dissipation induced by tidal torques (Sect. 4) in both 
short and long period dwarf novae (thus for both small 
and large accretion discs) using two characteristic sets of 
parameters taken from Ritter & Kolb (1998): those of SU 
UMa and SS Cyg (Table 0). For the short period sys- 
tems, although we take the parameters of SU UMa, we 
include neither the tidal instabilities (suggested by Osaki 
1989) to be responsible for superoutbursts in SU UMa), 
nor the combination of irradiation, mass transfer varia- 
tion and evaporation used by Hameury, Lasota & Warner 
(2000) to obtain long lasting outbursts. Therefore, we do 
not expect to reproduce superoutbursts. We also neglect 
irradiation of secondary star which in principle may have 
significant impact on some of our results. We will consider 
this effect in a future work. 





Long period 


Short Period 


Mi/M© 


1.2 


0.8 


M2/M0 


0.7 


0.15 


Porh (hr) 


6.6 


1.83 


rcirc (10^° cm) 


1.14 


0.91 


< rout > {W^° cm) 


5.4 


2.3 



3. Stream impact 

3.1. Description 

The gas stream from the secondary hits the accretion disc 
creating the observed bright spot. As the dynamical time 
is smaller than the thermal time (tdyn/^th = a < 1), the 
temperature is roughly constant around each annulus of 
the accretion disc. Thus, the effect of stream impact could 
be approximatively axisymmetric and we can use a ID 
approximation. This assumption is obviously very rough. 
Spruit & Rutten (1998) have shown that the stream im- 
pact region in WZ Sge is not at all axisymmetric (see also 
Marsh & Horne 1990 for IP Peg). Moreover, if we refer 
to the generally assumed values of a, idyn/^th is not that 
small in outburst; in addition, as the hot spot temperature 
Ths ^ Tdisc in the outer layers, the actual thermal time 
is much smaller than the unperturbed equilibrium value. 
Finally, we do not know exactly how the stream impact 
interacts with the disc and where precisely its energy de- 
posited in the disc; e.g. in the simulations of Armitage & 
Livio (1996, 1998) a significant amount of the stream ma- 
terial can ricochet off the disc edge and overflow toward 
smaller radii. 

It is therefore difficult to model properly the hot-spot 
contribution to the energy balance of the outer layers of 
an accretion disc. Here we consider that the stream im- 
pact heats an annulus Arhs of the disc with an 'efficiency' 
rji. The coefficient rji is supposed to represent the rather 
complex way the outer disc is heated by the stream. The 
energy per unit mass released in the shock between the 
stream and the disc is (V2 — Vi)^, where V2 is the kep- 
lerian velocity at the outer disc edge, IV2I — GAfi/(2rout), 
and Vi is the stream velocity at the position of impact. 
As the angle between Vi and V2 is large, and | Vi | is com- 
parable to IV2I, IV2 - Vil is of order of GMi/{2rout)- 

The heating rate Qi is therefore taken as: 



Qi{r) - V 



GM1M2 



1 



2rout 27rroutArh 



■ exp 



^out 



Arh, 



(1) 



where Mi is the primary's mass, M2 the mass transfer rate 
from the secondary and Tout the accretion disc outer ra- 
dius. This assumes that the difference between the stream 
and the Keplerian kinetic energy is released in a layer of 
width Arhs with an exponential attenuation. Note that 771 
can be larger than unity, as IV2 — Vi| can be larger than 
the Keplerian velocity. In the following, we take 771 = 1. 
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The energy equation then becomes: 

2(g+ + l/2Qi - Q- + J) UT^ldrVr 



dt 



fiCp r 



dr 
dTc 
dr 



(2) 



Here, Tc is the central temperature, Q'^ — (9/8)t^S]f2|^ and 
Q~ = crTJg the heating and cooling rate respectively, S 
the surface density, v the kinematic viscosity coefficient, 
J is a term that accounts for the radial energy flux carried 
either by viscous processes or by radiation (see Hameury 
et al. 1998), and Vr is the radial velocity. Cp is the specific 
heat at constant pressure, /i the mean molecular weight, 
and TZ the perfect gas constant. 

The cooling rate Q~(S, Tc, r) is obtained from interpo- 
lation in a grid obtained by solving the vertical structure of 



the disc (Hameury et al 1998). These structures are calcu- 



lated assuming a steady state, the departure from thermal 
equilibrium being accounted for by an effective viscosity 
parameter acs different from the actual viscosity parame- 
ter a. We assume here that the effect of the stream impact 
on the vertical structure can also be accounted for in this 
way. This would be perfectly correct if the local heating 
term were proportional to pressure, as is viscosity; if on 
the other hand it were localized at the disc surface, the 
problem would resemble that of disc illumination, and also 
in such a case, the use of an effective acs leads to qualita- 
tively correct results: Stehle & King (1999), who used this 
approximation, and Dubus et al. (1999), who calculated 
the exact vertical structure of irradiated discs, obtained 
similar modifications of the S curves by irradiation of the 
disc. 

3.2. Results 

Fig. H displays the E — Toff curves, calculated at two dif- 
ferent radii with and without the heating by the stream 
impact. These are the well known 'S-curves' in which 
three branches can be distinguished: the stable, low- and 
high-state branches where T^g increases with E, and the 
unstable middle branch where T^s decreases with E in- 
creasing. The latter branch is delimited by two critical 
points (Ei„in,T'min) and (E^ax, Tmax)- To compute these 
S-curves, a steady state was assumed so Toff is propor- 
tional to M^/^. As mentioned above, we take r]i — 1. 

When stream impact is taken into account one ob- 
serves at large radii significant deviations from the stan- 
dard case. It must, however, be stressed that Q\ is a non- 
local quantity proportional to Mi and thus depending on 
the outer boundary condition, whereas the usual S-curves 
are locally defined. Hence, the interpretations of these 
curves is not straightforward. For a fixed E, there is a fac- 
tor A log Toff ^ 0.5 between curves 2 and 2 . For a fixed 
Toff, S is reduced in the outer disc region when heating 
by the stream impact is included. It should therefore be 
easier, for a fixed mass transfer rate, to trigger outbursts 




Fig. 1. E— Teff curves for the parameters of SS Cyg . Curve 
1 is calculated at a radius 3.75 x 10^° cm that is out of the 
stream impact region. Curves 2 and 2 are calculated at 
a radius 5.2 10^° cm with and without the stream impact 
included in the model. For the hot branch Ofhot =0.2, for 
the cold one acoid = 0.04. The S-curve is obtained using 
Eq. (40) in Hameury et al. 1998 




Fig. 2. Emin (lower curves) and Emax (upper curves) for 
SS Cyg parameters. The dot-dashed lines show the results 
of the standard DIM, the solid lines have been obtained 
when the stream impact is taken into account, and the 
dashed lines when tidal effects are included. 



starting at the outer edge of the disc when heating by the 
stream impact is included. One also expects the disc to be 



less massive. Moreover, the Ej^ 



interval decreases 



(see Fig. 
too. 



so that the recurrence time should decrease 
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Fig. 3. Comparison of light curves with (solid line) and 
without the stream impact effect (dashed line). The pa- 
rameters we use are taken from table |l| and M2 = 2 x 10^^ 
g s~\ ?7i = 1 and Arhs = fout/lO. 



Fig. H shows light curves calculated without and with 
heating for a long period system ("SS Cyg"). One can 
see that the recurrence time is reduced in the model with 
the stream impact and that the quiescent luminosity is in- 
creased, as expected, since the stream impact heats up the 
outer parts of the disc which dominate the quiescent lu- 
minosity. During quiescence, however, the luminosity still 
increases in contrary to observations. Figures ^ and || show 
how the outburst type and the recurrence times depend 
on the mass-transfer rate. Including the additional heating 
by the stream impact implies that outside-in outbursts are 
triggered for lower mass transfer rates than in the stan- 
dard DIM. Also the upper critical rate Mc, above which 
the disc is in a steady hot state, is lower as predicted by 
Meyer & Meyer-Hofmeister (1983). 

For larger Arhs, the quiescent luminosity is smaller, 
because the additional amount of energy per unit surface 
decreases in the outer parts of the disc (note also that the 
integral of Qi over the disc surface depends weakly on Arhs 
because of the simplified form of Eq. |l| where terms of or- 
der of Arhs /''out have been neglected). For smaller 771, the 
recurrence time and the quiescent magnitude are higher. 
As short period systems have small discs, one might have 
expected that the stream impact would have stronger ef- 
fects in these systems, but the results are really similar as 
can be seen in Figs. I and |. However, in these calculations 
we assumed that the average disc outer radius is equal to 
the tidal radius. This is an overestimate of the disc's size, 
especially for the SU UMa stars. Realistically small discs 
will be discussed in Sect. ^. 

The main effects of the stream-impact heating seen on 
these figures are: (1) The critical accretion rate Mc (Smin) 
above which the disc is stable, is decreased by a factor 
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Fig. 4. Radius rinit at which the instability is triggered 
and quiescence time as a function of the mass transfer 
rate from the secondary. The parameters are those of SS 
Cyg. Symbols correspond to the following: squares for the 
standard DIM; circles when the stream impact is included; 
triangles when the stream impact and the tidal dissipation 
are both included. The high mass-transfer end of each se- 
quence of symbols corresponds approximately to the sta- 
bility limit. 



of more than 3 ; (2) the transition between inside-out and 
outside-in outburst ('A-B transition') occurs for a transfer 
rate Ma-b that is decreased by a factor < 2 ; (3) the 
recurrence time is almost unchanged. As expected (Osaki 
1996), the quiescence time is approximately constant when 
the outbursts are of the inside-out type, and decreases 
(very roughly as M2~^) ^'^^ outside-in outbursts, although 
for the smaller disc some deviations from this relation can 
be seen (see Sect. 6). (Note that we plot here the time in 
quiescence and not the recurrence time.) 

For systems with parameters close to that of SS Cyg, 
the inclusion of the stream impact heating lowers the crit- 
ical rate for outside-in outburst from Ma-b — 2.5 x 10"'^'' 
g s~^ to Ma^b — 1-7 X 10^^ g s^^ . The quiescence time 
is also shorter but only by a few days as long as M2 is not 
too close to Mc (see Fig. ||). 

For short period systems, the critical mass transfer is 
Mc ~ 5 X 10^^ g s~^, however, the transition between 
inside-out and outside-in outbursts still occurs for mass- 
transfer rates {Ma-^b — 2 x 10^^ g s~^) much higher that 
what is estimated for these systems (less than a few x 10"'^^ 



g s see e.g. Warner 1995). We will discuss this problem 
in Sect. O. 



4. Tidal dissipation 
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Fig. 5. Same as Fig. ^ for short period systems 



4-.1. Description 

Accretion discs in close binary systems are truncated by 
the tidal effect of the secondary star at a radius close 
to the Roche radius, the so-called tidal truncation ra- 
dius: angular momentum transported outwards by viscous 
torques is returned to the binary orbital motion by the 
tidal torques acting on the disc's outer edge. In the stan- 
dard DIM model, tidal torques can be included in the 
angular momentum equation in the following way (Smak 
1984; Hameury et al. 1998): 



9S _ 3 
'dt ~ ^rd^^^2^ 



1 

27rr 



tid 



(3) 



where the tidal torque per unit of area can be written as 



Ttid = cri/Silorb - 



(4) 



where a is the orbital separation of the binary and c a nu- 
merical constant that allows to adjust the tidal truncation 
radius < rout > that the disc would have had if it were 
steady (this gives an average value of the outer radius). 
The r" term used in this formula determines the disc frac- 
tion that is significantly affected by the tidal-torque dissi- 
pation (it corresponds to radii where the Smin and Smax 
curves are decreasing with radius (Fig. ||)). We use n = 5 
which is quite an uncertain value, and we discuss in Sect. 
6 the effects of changing it. 

As pointed out by Smak (1984) and discussed 
in detail by Ichikawa & Osaki (1994), tidal torques 
(as viscous torques) induce a viscous dissipation 
(Papaloizou & Pringle 1977): 



Qtid = (^^K — f^orb)7t 



tid 



(5) 



This tidal dissipation is usually neglected as it is small 
compared to except near the outer radius. During 



outbursts, the disc radius varies significantly, so does the 
tidal-torque dissipation which can be large enough to mod- 
ify the outburst behaviour. We therefore include this term 
in the energy equation: 



dt 



2{Q+ + l/2Qtid + l/2Qi -Q- + J) 



7^^c 1 drv, 
fiCp r dr 



a7 



4-2. Standard DIM plus tidal dissipation 

We first discuss the tidal effect alone, without including 
the heating by the stream impact. We find that type A 
outbursts are obtained for lower mass transfer rates than 
in the standard case; Ma-b is only slightly larger than 
in the case where the stream impact is included. For long 
period systems, we found Ma-b — 1-8 x lG^''g s~^. The 
critical mass transfer rate Mc is also slightly larger than 
in the model with stream impact. 

The major difference with previous calculations is the 
presence of a critical mass transfer rate AfgL , for which one 
obtains a sequence of alternating short (narrow) and long 
(wide) outbursts with roughly the same amplitudes (see 
Fig. |g) . At low mass transfer rates only short outbursts are 
present whereas for mass transfer rates higher than Msl 
all outbursts are long. In a small range of mass transfer 
rates outbursts of both durations are present (see Fig. ^ . 
We return to this point in the next section. 

Msl is comparable to A/a-b; for the parameters used 
in Fig. ^, AfsL < AfA-B, and short outbursts are all of the 
inside-out type. They are similar to those produced in the 
standard model, except that the quiescent flux is slightly 
larger as a result of viscous dissipation in the outer disc (as 
in the case of stream impact). As for the stream-impact 
heating case, the outburst peak luminosity is very slightly 
lower than in the standard model for the same mass trans- 
fer rate. The recurrence time scales are not affected by 
tidal dissipation when only short outbursts are produced. 
Long outb ursts are slightly brighter than the shor t ones as 
observed ( Oppenhcimer, Kenyon fc Mattel 1998 ) , and for 
this set of parameters they can be either of the inside-out 
or outside-in type because the critical mass-transfer rate 
above which they occur AfsL is less than A^a-b- 

Long outbursts appear because near the outburst max- 
imum the outer radius increases, so does tidal dissipa- 
tion. Tidal dissipation lowers the value of Smin so that 
it takes longer to reach this critical density at which a 
cooling-front starts to propagate and shuts off the out- 
burst. Outbursts therefore last longer than in the stan- 
dard case and have the characteristic "flat top" shape. The 
amplitude of the outer-disc variations increases with the 
mass-transfer rate. At low mass-transfer rate these varia- 
tions are too small to produce an observable effect on the 
outburst properties. The critical rate Msl corresponds to 
the case when the disc expansion begins to affect the out- 
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250 
time (days) 

Fig. 6. Light curve for the DIM plus tidal dissipation 
(without stream impact) for M2 = 1.5 x lO^^g s^^ and 
binary parameters of SS Cyg. The Both short (narrow) 
and long (wide) outbursts are of inside-out type. 



burst duration. Short outburst follow long ones because 
the tidal dissipation lowers both the critical and the ac- 
tual post-outburst surface-density so that during the next 
outburst the time to empty the disc is shorter. Sequences 
of alternating short and long outbursts are present for 
rather narrow interval of mass-transfer rates. At higher 
mass-transfer rates outbursts are of the outside-in type 
and therefore long (see Hameury et al. 1998). 

Compared to the standard case, the recurrence time is 
increased when long outbursts are present, but the dura- 
tion of quiescence is unchanged: the duration of the out- 
burst is the only modified quantity. 

When stream impact heating and the tidal dissipative 
process are both included in the model, their effects cu- 
mulate: MsL, Ma-b and Mc are smaller by about 15 % 
than what one gets when including only one effect (see 
Fig. I for Ma-b and Mc). 

5. The nature of long and short outbursts 

We have therefore a natural explanation of the coexistence 
of short and long outbursts in long period (> 3 hr) sys- 
tems: for which the mass transfer is large enough for the 
additional heating to produce this effect; in short period 
system, M2 corresponds to short outbursts only. Heat- 
ing by stream impact and tidal-torque dissipation account 
also for the occurrence of outside-in outbursts in SS Cyg, 
provided that the mass transfer rate happens to be large 
enough (approximatively two to three times larger than 
average) . 

In many cases, a bimodal di stribution of the outbursts 
is observed ( ^an Paradijs 198^ ). In the case of SS Cyg, 
Cannizzo & Mattel (1992) identified short outbursts with 



a duration between 5 and 10 days, and long ones with 
a duration between 13 and 17 days. Cannizzo (1993) ob- 
tained in his model a sequence of alternating short and 
long outbursts, but this was due to the use of a fixed 
outer radius which is an incorrect ('brick wall') boundary 
condition ( Hameury et al 1998| , ^mak 19*9^ ). In any case 
his sequence looks different from the observed one and all 
his outbursts are inside-out (type B). 

When tidal dissipation is included, one obtains, for 
a SS Cyg-like system, an alternating sequence of nar- 
row and wide outbursts for a mass transfer rate M2 ~ 
MsL ^ 1.5 X 10^''g (the value used in Fig. ||). As 
the short /long sequence of outbursts is not uncommon, a 
small mass transfer rate variation around MgL is the most 
likely cause for the occurrence of both narrow and wide 
outbursts. Moreover, if one believes that both A and B 
outburst types occur in SS Cyg (or another long period 
system), M2 must also vary around Ma^b- 

The value of A/gL depends on the binary parameters 
and for a given orbital period one obtains a small range 
of mass transfer rates for which both short and long out- 
bursts occur. Fig. 1^ shows this range as a function of the 
orbital period for a standard dwarf nova with a primary 
mass Ml = 0.6Mq and a main-sequence secondary mass 
M2 = 0.11 X Pill- M0. Note that now the primary's mass 
is different from the values used previously. Mi = 0.6Mq 
corresponds to the average white-dwarf mass in CVs. As 
mentioned earlier, < rout > is taken from a table in 
Paczynski (1977). The lower boundary of this range is Msl 
which varies by a factor of about 4 from a 2 hour period to 
a 7 hour period system. For short period systems, the mass 
transfer rate required to obtain both narrow and wide out- 
burst is larger than 2 x lO'^^g s~^. This is much larger than 
the expected secular values of a few times lO^^g s~^ and 
can explain why narrow and wide outburst alternation 
is not observed in short period dwarf novae. Superout- 
bursts in SU UMa's, which are different from long out- 
bursts (Warner 1995; Smak 2000), occur for mass transfer 
rates lower than Msl- On the other hand there exist a sub- 
class of the SU UMa stars, the ER UMa systems, which 
have very short supercycles, and are in outburst most of 
the time. They are presumed to have mass transfer rates 
as large as lO^^g s~^, of order of or above the value of 
Msl corresponding to their orbital periods, and should 
therefore show the alternation of long and short outbursts 
described here. This effect could add to the enhancement 
of the mass-transfer rate resulting from the illumination of 
the secondary that also cause sequences of long and short 
outbursts in short period systems (Hameury et al. 2000). 

Fig. shows also A/a-b; in this case we find that AfgL 
is larger than Ma-b, which means that all long outbursts 
are now of the outside-in type. However, the resulting light 
curves in which an alternance of short and long outbursts 
is found do not differ much from the one shown in Fig. ^ 
except that the rise times are shorter. 



8 



Buat-Menard et al.: The nature of dwarf nova outbursts 



We propose that stochastic variations of the mass- 
transfer rate for the secondary cause the sequence of al- 
ternating short and long, as well as type A and type B 
outbursts. For the parameters of SS Cyg, we need mass- 
transfer rate variations in the range ~ 1 — 2 x 10^^ g s~^ 
to obtain the four types in the same light curve; this range 
depends on the system parameters, and in particular on 
the orbital period (see Fig. and on the mass of the pri- 
mary. Because outburst duration depends mainly on the 
disc's size (Smak 1999b) for short outbursts whose dura- 
tion is independent of the mass transfer rate we obviously 
reproduce the correlation (van Paradijs 1983) between the 
outburst width and the orbital period. It is much more dif- 
ficult to reproduce the observed correlation for long out- 
bursts, since their duration is quite sensitive to M2, whose 
variation with the orbital period is not well known. One 
must, however, be cautious when interpreting this correla- 
tion, as (i) the statistics is small, and (ii) this correlation 
may result in part from the very definition of "long" out- 
bursts: they must last longer than short outbursts, whose 
duration increases with the orbital period simply because 
of the increase of the disc's size (Smak 1999b). In the 
model by Smak (1999a), long outbursts have to be inter- 
rupted by M2 returning to a low value since otherwise the 
disc would remain steady. This requires the action of an 
external physical mechanism with an unspecified free pa- 
rameter. If for long outbursts the correlation is real, in 
Smak's model the mechanism terminating long outbursts 
would have to depend on the orbital period. 

Since Msl is much larger than the mass-transfer rates 
expected in systems below the period gap, we do not ex- 
pect long outbursts for periods shorter than 2 hours. This 
agrees with observations. Superoutbursts observed in SU 
UMa systems must have a different explanation; they are 
either due to a tidal instability, as proposed by Osaki 
(1996), and/or to the irradiation of the secondary, which 
is quite efficient in these short period systems (Hameury 
et al. 2000). 

We still cannot, however, account for the presence of 
outside-in outbursts in short period systems, whereas nor- 
mal outbursts in SU UMa systems are generally believed 
to be of this type (Osaki 1996). There are two possi- 
ble ways out of this difhculty: first, as discussed in the 
Sect. 6, SU UMa normal outbursts occur in truncated 
discs and could be of the inside-out type; second, the disc 
outer radius could be much smaller than what the stan- 
dard DIM predict (Smak 2000); indeed, observations of 
eclipsing dwarf novae by Harrop-AUin and Warner (1996) 
tend to substantiate this idea. The second option is also 
favoured by the 'direct' observations of outside-in out- 
bursts in two eclipsing SU UMa systems (see Sect. |6.2| 
for references). 




1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 

Porb (hour) 

Fig. 7. Range of mass transfer rates for which both long 
and short outbursts occur and critical rate, as a function 
of the orbital period. Here, Mi = 0.6 M© and M2 = 0.11 x 



Porb Mq. The critical rate = 9.5x lO^^M- 



(see 

Hameury et al 1998). For low M, only short outbursts are 



-0.899.68 

'10 



observed; above the upper curve, all outbursts are long. 
The upper solid line shows the critical mass transfer rate 
above which the disc is stable. 



6. Inside-out versus outside-in outbursts 

Outbursts of SS Cyg are believed to be of both A and 



B types (5mak 1984, Cannizzo, Wheeler fc Polidan 1986 



Mauche 1996, 3mak 1998). Thus, the mass transfer rate 
from the secondary probably varies around the A- 
B transition. In models, outside-in outbursts are usu- 
ally obtained only for mass transfer rates large com- 
pared to those expected from binary evolution models 



(Baraffe & Kolb 1999). Moreover, disc model computa- 
tions with such values of M2 predict recurrence times 
that are much shorter than observed. Ichikawa & Osaki 
(1992) managed to get rid of this problem by using an 
r-dependent a-prescription specially designed to suppress 
type B outbursts. Their viscous diffusion time is then al- 
most independent of the radius and inside-out outbursts 
occur only at very small mass transfer rates. 

As mentioned earlier, including the stream- 
impact heating lowers Ma-b without changing the 
a-prescription. We obtain Ma-b = 1-7 x 10^^ g s""'^, 1.4 
times smaller than in the standard case; this is reasonably 
close to estimates of SS Cyg average mass transfer rate 
M2 = 6 X lO^'^g s~^ (e.g. [Patterson 1984|). An increase 



of the mass transfer rate by a factor 2 - 3 is therefore 
sufhcient to provoke outside-in outbursts. On the other 
hand, for short period systems, the transition between 
inside-out and outside-in outburst still occurs for much 
too large mass transfer rates (2 x 10^^ g s~^), whereas 
estimates of mass transfer in thes e systems do not exce ed 
a few xlO^s g s-i (Warner 1995; [Baraffe fc Kolb 1999|) . 
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6.1. Predictions of the DIM 

Outside-in outbursts occur if the time scale for matter 
accumulation at the outer disc edge iaccum is shorter than 
the viscous drift time scale tdrift- These times are (Osaki 
1996): 



M2 



t. 



drift 



(7) 



(8) 



where v = 2/3a{TZ/ fi)T/Q is the kinematic viscosity, and 
S 0.05 is a numerical correction factor. Using the an- 
alytical fits of Smax given by Hameury at al. (1998), the 
condition for outside-in outbursts is: 



M > 4.35 X 10^^ 



T 



4000 K 



^-0.88^2^4^3-1 



(9) 



where rio is the radius measured in 10^" cm and where 
we have assumed TZ/fJ, — 5 x 10''. This condition can be 
compared with our numerical results. In the case without 
the additional effects, our critical rate is 2.5 x 10^^ g s~^ 
for long period systems, whereas Eq. gives 1.8 x 10^'' 
g s~^; for short period systems, these two numbers are 

2.5 X 10^^ and 2.8 x 10^^ g s~^ respectively. It must be 
noted that the quality of the agreement should not be 
a surprise, since the quantity S was introduced by Osaki 
(1996) precisely for obtaining good fits. 

When heating of the outer parts of the disc is included 
Eq. (^ is no longer valid. Osaki (1996) assumed that in 
quiescence the surface density is parallel to the critical 
density profiles, that S 2I]min, and he used the stan- 
dard critical densities. Things arc different when heating 
of the outer disc is included. Near the disc outer rim both 
Emin and Emax are considerably reduced as is the distance 
between them. These effects reduce the value of Ma-b- 

This reduction is not sufficient to obtain desired values 
of the mass-transfer rate, i.e. to ~ 3 x 10^^ g s~^, which 
is typical of the mass transfer rate in SU UMa systems. 
The only hope is in the strong dependence of Ma-b on 
the outer disc radius. The mean radius we have assumed 
for short period systems (2.3 x 10^" cm) is probably too 
large for most SU UMa's since it exceeds in many cases 
the 3:1 resonance radius, which is of order of 0.46 a, i.e. 
2 X 10^*^ cm for a total mass of 0.7 M0 and a period of 

1.6 hr characteristic of SU UMa stars. The 3:1 resonance 
radius should be reached only during superoutbursts and 
should be much smaller (0.3 a, i.e. 1.3 x 10^° cm) during 
the first few outbursts of a supercycle . For such small radii 
Ma-b would then be reduced by a factor 4 - 5, i.e. to a 
value slightly larger than the estimated mass transfer rate 
in SU UMa systems. 

In Fig. 1^ we show the quiescence time and the radius 
at which the outburst starts versus mass transfer rate, 
for a short period system with a mean outer disc radius 
~ 1.3 X 10^° cm. This model corresponds to a relaxed 



disc, i.e. a disc whose mass is constant over an outburst 
cycle. Such small discs are expected to exist only during 
the first half of a supercycle and later should grow in size. 
We have assumed a large coefficient c in Eq. ^ in order 
to obtain a small disc, but we have not included heating 
by tidal torques that would have resulted in much too 
large a dissipation: the dissipation of tidal torques must 
be negligible in discs much smaller than the truncation 
radius. As can be seen, outside-in outbursts occur for mass 
transfer rates as low as 5 — 6 x 10"'^^ g s""'^, reasonably close 
to the expected value in these systems. 

It is worth noting that the recurrence time depends 
on the mass transfer rate, even in the case of inside-out 
outbursts. A comparison of Figs. ^, ^ and || shows that 
this effect is related to the size of the disc. This is due to 
the fact that a local dimensional analysis becomes invalid 
when the radius of the outer edge of the disc is not much 
larger than the inner radius (for example, in the case of 
steady discs, the inner boundary condition introduces a 
factor [1 — inn/ry^'^])- Osaki (1996) asserts that the cor- 
relation tg oc found between the supercycle recurrence 
time is and the recurrence time of normal outbursts t^ 
in SU UMa systems results from Eq.(|^) and the relation 
ts oc which results from the assumption that super- 

outburst occur when sufficient mass is accumulated (Os- 
aki 1996; Hameury et al. 2000). Clearly, this cannot be the 
case if one uses the standard disc model in which a is not 
forced to adopt a functional dependence that would give 
the desired effect but is just kept constant. In any case, 
the mass transfer rate cannot be assumed to be the only 
variable on which recurrence times depend. For example 
Menou (2000) showed that both the normal outburst and 
super-outburst recurrence times strongly depend on the 
secondary to primary mass ratio. 

Small discs are also expected if one reduces the index n 
in Eq. (||), since the smaller n, the larger the amplitude of 
disc radius variations during outbursts; for a given maxi- 
mum disc radius (some fraction of the Roche radius), the 
disc radius during quiescence is therefore expected to be 
smaller for smaller n. 



6.2. Observations 

Campbell (1934) identified four distinct types of outburst 
in SS Cyg light curves which have a large variety of char- 
acteristic rise times. Since then, several techniques have 
been used to investigate the physical outburst properties: 
multi-wavelength photometry, spectrometry, eclipse map- 
ping and model fitting. The first disc instability models 
have predicted the outburst types A (outside in) and B 
(inside-out) ( ^mak 1984 ), and showed that the first cate- 
gory rises earlier at optical wavelengths and have a more 
asymmetric shape. As a result, asymmetric outbursts are 
usually supposed to be of type A. However asymmetric 
outbursts can also be of the inside-out type. For instance. 
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Fig. 8. Radius i?init at which the instabihty is triggered 
and quiescent time as a function of the mass transfer rate 
from the secondary. The parameters are those of short 
period systems with < rout >~ 1-3 x 10^° cm. Squares 
correspond to the standard DIM, and circles to the case 
where the stream impact has been inchided. 



in Fig. ^, long outbursts are not symmetric but are of type 
B. 

The UV delay is the observed lag between the rise in 
optical and the rise in UV. Its duration is believed to be a 
signature of the outburst starting region. In the naive pic- 
ture since the hottest parts of the disc associated with the 
shortest wavelengths are close to the white dwarf, there 
should be almost no UV delay if the outburst starts in 
the inner part of the accretion disc, whereas one should 
expect an important UV delay for outside-in outbursts. 
However, as discussed by Smak (1998), things are more 
complicated than what the naive picture would suggest, 
and long delays can also exist for inside-out outbursts if 
the outer disc radius expansion during outburst and the 
correct size of the disc are taken into account. For a given 
binary system the UV delay will always be longer for an 
outside-in outburst than for an inside-out one, but the 
length of the delay is not sufficient by itself to determine 
the outburst type. 

Smak (1998) proposed that the ratio of the UV delay 
and the outburst duration is a good indicator of the out- 
burst type. However, the existence of both long and short 
outbursts at a given orbital period makes this test very 
uncertain. Indeed Fig. |^ shows a probing counter-example 
where both short and long outbursts are of type B and for 
which Smak's ratio would be quite different. 

In addition Hameury et al. (1999), using SS Cyg pa- 
rameters, reproduced the observed ( [Mauche 1996| ) asym- 
metric outburst and 1-day EUV delays for an inside-out 
outburst starting to propagate at the inner edge of a trun- 



cated disc. This obviously raises the question of the real 
relationship between the outburst type and the UV delay. 

It has been known for some time that eclipse profiles 
allow to determine the propagation direction of an out- 
burst (Smak 1971 for U Gem; Vogt 1984 and Rutten, van 
Paradijs &Tinbergen 1992 for OY Car). For example in 
the case of an accretion disc whose luminosity is much 
larger than both the white dwarf and the hot spot lumi- 
nosity, the eclipse profile is almost symmetric. This is the 
case of HT Cas, a SU UMa system with a period of 106 



minutes (loannou et al. (1999)). The authors found that 
during the rise of the burst, the eclipse is shallow and its 
width is large. This means that most of the flux origi- 
nates from the outer parts of the disc. They observed that 
the width of the eclipse decreases after the outburst peak, 
which is in perfect agreement with the predicted decrease 
of the disc radius during decline, therefore proving the 
outside-in nature of the outburst. Webb et al. (1999) and 
Baptista, Catalan & Costa (2000) used the same technique 
to bring evidences of the inside-out nature of an outburst 
of IP Peg and EX Dra respectively. 

At present, the only way to be sure of the outburst type 
is to use, whenever possible, eclipse profiles. As eclipses are 
not present in all systems it would be worthwhile to test 
Smak's assumption on the UV delay in eclipsing dwarf 
novae. 

7. Conclusions 

We have shown that the inclusion of additional heating 
effects such as the stream impact and tidal dissipation 
might solve several problems of the standard DIM. 

The most important result is that when the tidal 
torque dissipation in the outer disc regions is in- 
cluded, one obtains for a certain range of mass- 
transfer rates an alternation of short and long out- 
bursts which is observed in long orbital period sys- 
tem such as SS Cyg (Cannizzo & Mattel 1992) or Z 
Cam (Oppenheimer, Kenyon & Mattel 1998). The long 



outbursts are slightly brighter than short outbursts as ob- 
served. Moderate fluctuations of the mass-transfer rate 
would extend the range for which such alternating out- 
bursts are present. For systems with mass transfer rates 
close to the stability limit, such fluctuations will produce 
standstills as shown in Buat-Menard et al. (2000). This 
solution of the long/short outburst problem (Smak 2000) 
is different from the one suggested by Smak (1999a). In 
our model mass-transfer fluctuations play only an auxil- 
iary role and long outbursts do not result from fluctuations 
bringing the disc to a steady state. 

We showed that additional heating effects make the 
outside-in outburst possible at low mass transfer rates 
without changing the a-prescription. Our results on mass 
transfer rates and outburst types are in good agreement 
with most estimates in the case of long period systems. 
In the case of short period systems of the SU UMa type. 
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for outside-in outburst to occur the disc would have to be 
smaller than usually assumed. For small discs the recur- 
rence times of both the inside-out and outside-in outbursts 
is M2-dependent. It must be stressed out that there are 
still many uncertainties on the observational determina- 
tion of outburst types. For instance, the UV delay is not 
necessarily a good criterion for determining the outburst 
type; the eclipse profile during outburst is a much better 
criterion, but unfortunately restricted to a small subset of 
dwarf novae. 

There are still several points that need to be clarified. 
First, the stream impact effect geometry is not well known 
and is considered here in a ID approximation (as for the 
tidal dissipation). We do not know how far and how effi- 
ciently the energy propagates in the disc. The expression 
of the tidal torque can also be subject to discussion. More- 
over, we did not include irradiation here, and it would be 
interesting to study the combined effects of all phenom- 
ena. Finally, as raised above, the existing discrimination 
criteria on the outburst types A and B have to be tested 
and improved. 
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